N95- 22992 



*10 


COVARIANT HARMONIC OSCILLATORS AND 
COUPLED HARMONIC OSCILLATORS 



) 




D. Han 

National Aeronautics and Space Administration, Goddard Space Flight Center, 
Code 910.1, Greenbelt, Maryland 20771 

Y. S. Kim 

Department of Physics, University of Maryland, College Park, Maryland 20742 

Marilyn E. Noz 

Department of Radiology, New York University, New York, New York 10016 


Abstract 

It is shown that the system of two coupled harmonic oscillators shares the basic sym- 
metry properties with the covariant harmonic oscillator formalism which provides a concise 
description of the basic features of relativistic hadronic features observed in high-energy lab- 
oratories. It is shown also that the coupled oscillator system has the SL(4,r) symmetry in 
classical mechanics, while the present formulation of quantum mechanics can accommodate 
only the Sp(4,r) portion of the SL(4,r) symmetry. The possible role of the SL(4,rj symmetry 
in quantum mechanics is discussed. 


1 Introduction 

The covariant harmonic oscillator formalism developed by the present authors has been shown 
to be effective in explaining the basic phenomenological features of relativistic extended hadrons 
observed in high-energy laboratories. In particular, the formalism shows that the quark model and 
Feynman’s parton picture are two different manifestations of one relativistic entity. In addition, 
the formalism constitutes a representation of Wigner’s little group for a massive particle with 

internal space-time structure [1]. . . . , 

Since the classical mechanics of two coupled harmonic oscillators is discussed in Goldstein s 
text book [2], there is a tendency to believe that this oscillator problem is completely understood 
and that nothing new can be learned from it. We disagree. In this paper, we show that this 
coupled oscillator system can serve as an analog computer for the above-mentioned covariant 
oscillator formalism. 

From the mathematical point of view, the standard approach is to construct a suitable repre- 
sentation of the symmetry group after writing down its generators. The first symmetry group in 
the present case is 5p(4,r) with ten generators [3, 4, 5]. The second symmetry group is 51(4, r) 
which contains a number of 5p(4)-like subgroups. In constructing these groups, we shall note that 
each oscillator has its own 5p( 2) symmetry, and that the coupling of the two oscillator also has a 
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Sp(2)-like symmetry. It was pointed out that these three Sp( 2) groups can be combined into one 
Sp(4) group. 

Since Sp(4) is locally isomorphic to the deSitter group 0(3,2), it can explain the Lorentz 
transformation properties, particularly that of the covariant harmonic oscillator formalism. In 
this paper, we concentrate on the issue of a lack of information on one oscillator affecting the 
uncertainty and the entropy of the other oscillator. 


2 Covariant Harmonic Oscillators 


The covariant harmonic oscillator formalism has been discussed exhaustively in the literature, and 
it is not necessary to give another full-fledged treatment in the present paper. Instead, we shall 
concentrate on the issue of entropy in this paper. The entropy is a measure of our ignorance and 
is computed from the density matrix [6, 7]. The density matrix is needed when the experimental 
procedure does not analyze all relevant variables to the maximum extent consistent with quantum 
mechanics. The purpose of the present note is to discuss a concrete example of the entropy arising 
from our ignorance in relativistic quantum mechanics. 

Let us consider a bound state of two particles. For convenience, we shall call the bound state 
the hadron, and call its constituents quarks. Then there is a Bohr-like radius measuring the space- 
like separation between thd quarks. There is also a time-like separation between |;he quarks, and 
this variable becomes mixed with the longitudinal spatial separation as the hadron moves with a 
relativistic speed. 

However, there are at present no quantum measurement theories to deal with the above- 
mentioned time-like separation. We shall study in the present paper how this ignorance is trans- 
lated into the entropy. Within the framework of the covariant harmonic oscillator formalism [1], 
it will be shown that the entropy increases as the hadron gains its speed. The entropy defined in 
this way is a more fundamental quantity than the hadronic temperature [4], It is independent of 
the question of whether the temperature can be defined [8]. 

Let us consider a hadron consisting of two quarks. If the space-time positions of two quarks 
are specified by x a and x b respectively, the system can be described by the variables [9] 

X = (x a + x b )/2, x = (x a - x b )/2y/2. (1) 

The four-vector X specifies where the hadron is located in space and time, while the variable x 
measures the space- time separation between the quarks. In the convention of Feynman et al [9], 
the internal motion of the quarks bound by a harmonic oscillator potential of unit strength can 
be described by the Lorentz-invariant equation 

if 5 2 1 

2 H ~ 1 = A ^ X >‘ (2) 

We use here the space-favored metric: x ** = ( x,y,z,t ). 

It is possible to construct a representation of the Poincare group from the solutions of the 
above differential equation [1]. If the hadron is at rest, the solution should take the form 


/ 1 \ V4 

il>(x,y,z,t) = ij>(x,y,z) [ — ) exp 


(-*72) , 


(3) 
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where ^(x, y, z ) is the wave function for the three-dimensional oscillator with appropriate angular 
momentum quantum numbers. There are no excitations along the t direction. Indeed, the above 
wave function constitutes a representation of Wigner’s 0(3)-like little group for a massive particle 

[i]. 

Since the three-dimensional oscillator differential equation is separable in both spherical and 
Cartesian coordinate systems, tj>(x y y,z) consists of Hermite polynomials of z,y, and z. If the 
Lorentz boost is made along the z direction, the x and y coordinates are not affected, and can be 
dropped from the wave function. The wave function of interest can be written as 

tl> n {z, t) = ( J ) 1/4 exp ( -t 2 / 2 ) t/>„(2), (4) 


/ 1 \ 1 / 2 

r[z) = (^jF) H n ( z ) exp(— z 2 /2), (5) 

where ip n {z) is for the n-th* excited oscillator state. The full wave function is 

The subscript 0 means that the wave function is for the hadron at rest. The above expression is 
not Lorentz-invariant, and its localization undergoes a Lorentz squeeze as the hadron moves along 
the z direction [1]. 

It is convenient to use the light- cone variables to describe Lorentz boosts. The light-cone 
coordinate variables are 

u = (z + t)/\/ 2, v = (z — t)/\/ 2. (7) 

In terms of these variables, the Lorentz boost along the z direction, 


/ cosh 7} sinh rj 
\ sinh 7/ cosh 


;)(:)■ 


(8) 


takes the simple form 

u — e v u , v — (9) 

where 77 is the boost parameter and is tanh _1 (u/c). The wave function of Eq.(6) can be written 

1 x 1 / 2 


as 


«(*.<) = GsH?) «„((« + »)/V2)exp{-i(u 2 + l> ! )}. 


( 10 ) 


If the system is boosted, the wave function becomes 


<(M) = (^)' /! H » + f2V )} • (>*) 

As was discussed in the literature for several different purposes, this wave function can be expanded 

35 ^ 1/2 

V^(z,t) = (l/coshr?r +1 X; ~ ) (tanh r]) k ip n+k {z)xp n (t). (12) 
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In both Eqs. (10) and (11), the localization property of the wave function in the uv plane 
is determined by the Gaussian factor, and it is sufficient to study the ground state only for the 
essential feature of the boundary condition. Eq.(10) and Eq.(ll) then respectively become 

$?(M)=(^) exp j-^(u 2 + w 2 )J. ( 13 ) 

If the system is boosted, the wave function becomes 

exp j-i(e -2 V + e 2 V)J. ^ (14) 

We note here that the transition from Eq.(13) to Eq.(14) is a squeeze trasnformation. The wave 
function of Eq.(13) is distributed within a circular region in the uv plane, and thus in the zt plane. 
On the other hand, the wave function of Eq.(14) is distributed in an elliptic region. This ellipse 
is a “squeezed” circle with the same area as the circle. This Lorentz-squeezed wave function can 
be expaned as 

tMM) = ^^Xl(tanh r}) k ^ k (z)^ k (t). (15) 

From this wave function, we can construct the pure-state density matrix 


Pv( z > z \ 0 = i’v (*> OiM*'* t '), 

which satisfies the condition p 2 = p\ 

P,{z,t-z\t') = J p v (z 1 t-,z",t")p r) (z",t"-,z',t')dz"dt". 


(16) 


( 17 ) 


However, there are at present no measurement theories which accommodate the time-separation 
variable t. Thus, we can take the trace of the p matrix with respect to the t variable. Then the 
resulting density matrix is 


Pr,{ z , z ') = J iPniZy l ) {^v( z \ t)}*dt 

- (dM)V tanh,) “ w ^ M - 


(18) 


The trace of this density matrix is one, but the trace of p 2 is less than one, as 

Tr {p 2 ) = J Pv{ z , z ')pr,( z ', z ) dz ' dz (19) 

= (dh^) ? (Unh ’’>“ 

which is less than one. This is due to the fact that we do not know how to deal with the time-like 
separation in the present formulation of quantum mechanics. Our knowledge is less than complete. 
The standard way to measure this ignorance is to calculate the entropy defined as [6, 7] 

S = -Tr (p\n(p)) . 
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If we pretend to know the distribution along the time- like direction and use the pure-state density 
matrix given in Eq.(16), then the entropy is zero. However, if we do not know how to deal with 
the distribution along t , then we should use the density matrix of Eq.(18) to calculate the entropy, 
and the result is 

S 4= 2 {(cosh T]) 2 ln(cosh tj) — (sinh q) 2 ln(sinh 7) j . (20) 

In terms of the velocity v of the hadron, 


5 = — ln[l - (v/c) 2 ] - 


(v/c) 2 ln(u/c) 2 
1 - (v/c) 2 


( 21 ) 


We can also calculate the density matrix using the Gaussian form of the wave function given 
in Eq.(17), and the result is 





z + z') 2 / cosh 2t] + ( 


z — z') 2 cosh 2t/]| , 


( 22 ) 


This expression also leads to the entropy given in Eq.(20) 
The diagonal elements of the above density matrix is 


/ 1 y /2 

' (2 ’ z)= U^) ex PHVco S h2,) 


(23) 


The width of the distribution becomes (cosh r?) 1/2 , and becomes wide-spread as the hadronic speed 
increases. Likewise, the momentum distribution becomes wide-spread. This simultaneous increase 
in the momentum and position distribution widths is called the parton phenomenon in high-energy 
physics. The position- momentum uncertainty becomes coshf/. This increase in uncertainty is due 
to our ignorance about the physical but unmeasurable time-separation variable. 

The use of an unmeasurable variable as a “shadow” coordinate is not new in physics and is of 
current interest [10, 11, 12, 13]. Feynman's book on statistical mechanics contains the following 
paragraph [14]. 

When we solve a quantum- mechanical problem, what we really do is divide the universe into 
two parts - the system in which we are interested and the rest of the universe. We then usually 
act as if the system in which we are interested comprised the entire universe. To motivate the use 
of density matrices, let us see what happens when we include the part of the universe outside the 
system. 

In the present paper, we have identified Feynman’s rest of the universe as the time-separation 
coordinate in a relativistic two-body problem. Our ignorance about this coordinate leads to a 
density matrix for a non-pure state, and consequently to an increase of entropy. It is interesting 
to note that the density matrix of Eq.(22) becomes that of the harmonic oscillator in a thermal 
equilibrium state if (tanh rj) 2 is identified as the Boltzmann factor [15]. 

We have thus far studied the properties of covariant harmonic oscillators where the longitudinal 
and time-like coordinates undergo squeeze transformations. The word “squeeze is relatively new 
in physics. However, squeeze transformations are almost everywhere in physics. In the rest of this 
paper, we shall discuss the role of squeeze transformations in the system of two coupled harmonic 
oscillators. We shall see that the problem of covariant harmonic oscillators with two variables is 
the same as that of two coupled harmonic oscillators. 
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Linear Canonical and Non- Canonical Transformations 
in Classical Mechanics 


For a dynamical system consisting of two pairs of canonical variables X\,p\ and X2,/>2, we can 
introduce the four-dimensional coordinate system: 


(*?1, 7/2, f73,*?4) = (*l,X2,Pl,p 2 ). 

Then the transformation of the variables from 77, to £, is canonical if 

MJM = J, 

where 


(24) 


(25) 


and 


J = 



II 

6 , 


0 

0 

1 

0 \ 

0 

0 

0 

1 

-1 

0 

0 

0 

0 

-1 

0 

0 / 


For linear canonical transformations, we can work with the group of four-by-four real matrices 
satisfying the condition of Eq.(25). This group is called the four-dimensional symplectic group 
or Sp( 4). While there are many physical applications of this group, we are interested here in 
constructing the representations relevant to the study of two coupled harmonic oscillators. 

It is more convenient to discuss this group in terms of its generators G , defined as 


M = exp (—iaG), 


(26) 


where G represents a set of purely imaginary four-by-four matrices The symplectic condition of 
Eq.(25) dictates that G be symmetric and anticommute with J or be antisymmetric and commute 
with J. 

In terms of the Pauli spin matrices and the two-by-two identity matrix, we can construct the 
following four antisymmetric matrices which commute with J of Eq.(25). 


'-H 


0 < 7 j 

<Tx 0 


)• 


( 7 2 0 

0 cr 2 


). 


0 <r 3 

-cr 3 0 

The following six symmetric generators anticommute with J. 


)• 0 - 


(27) 


A, 


and 


(° *1, 

k 2 = l - ( 1 

° ) 

I< 3 = 

= --f 0 n), 


\<T 3 0 J 

2 VO 

-I) ’ 

2\a 1 0 /’ 


( -<?3 0 \ 
Vo aj' 

Ql = 2 (/ 

i)- 

g 3 = 

. i(* 1 ° \ 

2 V 0 -< 7 X J ' 

(28) 
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These generators satisfy the commutation relations: 

[ Ji,Jj ] = itijkJk, [«/,-, Aj] = ieijkhk, [At, A j] = — —if-xjkJki 

[A, Jo] — 0, [A",, Qj] = *<5,j J 0 , 

[A, Qj] = iti,kQk, [A’ n A] = [Q„ A] = —iKi- (29) 

The group of homogeneous linear transformations with this closed set of generators is called the 
symplectic group 5p(4). The J matrices are known to generate rotations while the K and Q 
matrices generate squeezes [4]. 

It is often more convenient to study the physics of four-dimensional phase space using the 
coordinate system 

( 6 , 6 > 6 > 6 ) = { x \iPu x iiPi) • ( 30 ) 

The transformation frorii (r) u r/ 2 , rj 4 ) is 


(31) 


fix \ 


/ 1 0 0 0 \ 



6 


0010 


V2 

6 , 


0100 


m 

uJ 


0 

0 

0 


\V 4 / 


and the J matrix becomes 

0 1 0 0\ 

-10 0 0 
J - 0 0 0 1 
0 0-10/ 

In this new coordinate system, the rotation generators take the form 

*-2(! '•)■ J -iC 


(32) 


Jrt 


2 \cr 2 0 
-1 ( a 2 


(a 2 0 \ j _ Z±( a ? 0 ") 

U — (T 2 ) J ° 2 \ 0 (72 / 


(33) 


The squeeze generators become 

*>- s ( 

q , A (-- 0 


o i 0 
0 — O , 


Ko = - 


(? :)■ o)- 


0 aj ' ^“2(0 <n)’ Q3 2(^3 o 3 )’ 


(34) 


In addition to the ten generators given in Eq.(33) and also in Eq.(34), we can consider the scale 
transformation in which both the position and momentum of the first coordinate are expanded 
and those of the second coordinate contracted. The Hamiltonian given in Eq.(46) suggests such a 
transformation, and the transformation can be generated by 


*4(J i)- 


(35) 
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This matrix generates scale transformations in phase space. The transformation leads to a radial 
expansion of the phase space of the first coordinate [16] and contracts the phase space of the 
second coordinate. What is the physical significance of this operation? As we discussed in Sec. 
7, the expansion of phase space leads to an increase in uncertainty and entropy. Mathematically 
speaking, the contraction of the second coordinate should cause a decrease in uncertainty and 
entropy. Can this happen? The answer is clearly No, because it will violate the uncertainty 
principle. This question, will be addressed in future publications. 

In the meantime, let us study what happens when the matrix So is introduced into the set of 
matrices given in Eq.(33) and Eq.(34). It commutes with Jo,J 3 ,K\,K 2 ,Q\, and Q 2 . However, its 
commutators with the rest of the matrices produce four more generators: 

T)' = I ( / o'). 

o). i s ”'^]= j(° 3 t). < 36 > 

If we take into account the above five generators in addition to the ten generators of 5'p(4), 
there are fifteen generators. They form the closed set of commutation relations for the the group 
SL(4, r). This SL(A,r) symmetry of the coupled oscillator system may have interesting physical 
implications. 


4 SL(4,r) Formulation of Two Coupled Oscillators 


Let us consider a system of two coupled harmonic oscillators. The Hamiltonian for this system is 

# =T (— p\ + — pl + A'x\ + B'x\ + C'x 1 x 2 } . (37) 

2 l m 1 m 2 J 

where 

A' >0, B' > 0, 4 MB' -C a > 0. (38) 

By making scale changes of x\ and :r 2 to {m\jm 2 y^x\ and (m 2 /mi) 1 ^ 4 a: 2 respectively, it is possible 
to make a canonical transformation of the above Hamiltonian to the form [17, 18] 

H = ^ {p\ + P 2 } + ^ {-4*1 + Bx\ + Cx!x 2 j , (39) 

with m = (m 1 m 2 ) 1 / 2 . We can decouple this Hamiltonian by making the coordinate transformation: 

(Vi \ _ f cos(a/2) - sin(a/2) \ / x, \ . . 

\ 3/2 / \ sin(a/2) cos(a/2) ) \x 2 )' 


Under this rotation, the kinetic energy portion of the Hamiltonian in Eq.(39) remains invariant. 
Thus we can achieve the decoupling by diagonalizing the potential energy. Indeed, the system 
becomes diagonal if the angle a becomes 


tan q = 


C 

B-A' 


(41) 
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This diagonalization procedure is well known 

We now introduce the new parameters K and tj defined as 


K = \fAB - C' 2 /4, exp (-2 77) = 


A + B t- /(A - Df + C 2 


v'4A/? - 


A = K ( e 2 ” cos 2 ^ + e 2 ’ ; sin 2 ^ ) , 

B = 1\ ( e 2 ” sin 2 


a _2 n 2 a \ 

-he ’’cos -), 


A = A' (e' 2 " - e 2 ”) s 
the Hamiltonian can be written as 


sin a. 


H = {«i + 9a} + T + e> ’ 


where y\ and ?/2 are defined in Eq (40), and 


q-i \ _ /cos(a/2) — sin((»/2) 


/ <7i \ _ / « 
V 92 / V si 


sin(o/2) cos(o/2) 7 VP2 


*7 2) W Pi 7 

/2) / V p2 ) 


(42) 


in addition to the rotation angle a. In terms of this new set of variables. 4, B ;u;d e take the 
form 


(43) 


(44; 


(45) 


This form will be our starting point. The above rotation together with that of Eq.(40 i is generated 

by Jo- 
lt we measure the coordinate variable in units of (m/7) 1 ' 4 , and use (m/\ ) i ' ! for the momentum 
variables, the Hamiltonian takes the form 


H = ^ (e-"9 2 + e”y 2 ) + V” (e^ 2 f e^y]) , 


(46) 


where w = yjKjm. If r/ = 0, the system becomes decoupled, and the Hamiltonian becomes 

»= £(»>?+*?) + f(rf+*0- (47) 

In Sec. 8, we will be dealing with the problem of what happens when no observations are made 
on the second coordinate. If the system is decoupled, as the above Hamiltonian indicates, the 
physics in the first coordinate is solely dictated by the Hamiltonian 

(?: + *?) ■ ( 48 ) 

It is important to note that the Hamiltonian of Eq.(47) cannot be obtained from Eq.(46) by 
canonical transformation. For this reason, the Hamiltonian of the form 

X'=j (e~\! + «*»?) + \ + e-jr J) (49) 

may play a useful role in our discussion. This Hamiltonian can be transformed into the decoupled 
form of Eq.(47) through a canonical transformation. 
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5 Quantum Mechanics of Coupled Oscillators 

It is remarkable that both the Hamiltonian H of Eq.(46) and H’ of Eq.(49) lead to the same 
Schrodinger wave function. If y x and y 2 are measured in units of ( mK ) 1 ^ 4 , the ground-state wave 
function for this oscillator system is 

M x u* 2 ) = ‘^exp -|- e _,, y 2 ) j. (50) 

The wave function is separable in the yi and y 2 variables. However, for the variables Xi and x 2 , 
the story is quite different. If we write this wave function in terms of X\ and x 2 , then 


V>(xi,z 2 ) = ^ ex P ^(*1 cos | ” x 2 sin ^) 2 


• a Or 2 1 1 

+e ^(xj sm — + x 2 cos — ) I >. 


If r] — 0 , this wave function becomes 





(51) 


(52) 


For other values of 77 , the wave function of Eq.(51) can be obtained from the above expression by 
a unitary transformation. 

) ( t 4 to1 nl 2 (®i (^"X )V*m 2 (*^ 2 )) (53) 

mj m 2 

where is the m th excited state wave function. The coefficients A mim2 (rj) satisfy the unitarity 

condition 

X) |4mim 2 (a,»?)| 2 = L ( 54 ) 

mi m 2 

It is possible to carry out a similar expansion in the case of excited states [1]. 

As for unitary transformations applicable to wave functions, let us go back the generators of 
canonical transformations in classical mechanics. As was stated before, they are also applicable 
to the Wigner phase-space distribution function. The canonical transformation of the Wigner 
function is translated into a unitary transformation of the Schrodinger wave function. There are 
therefore ten generators of unitary transformations applicable to Schrodinger wave functions [4, 3]. 

The Wigner phase-space picture is often more convenient for studying the problems of coupled 
harmonic oscillators. Unitary transformations in the Schrodinger picture can be achieved through 
canonical transformations in phase space. It has been known that canonical transformations are 
uncertainty-preserving transformations. They are also entropy-preserving transformations [5]. Are 
there then non-canonical transformations in quantum mechanics? 

In the present case of coupled harmonic oscillators, we assume that we are not able to measure 
the x 2 coordinate. It is often more convenient to use the Wigner phase-space distribution function 
to study the density matrix, especially when we want to study the uncertainty products in detail 
(18,14]. 
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For two coordinate variables, the Wigner function is defined as [18] 

W{x l ,x 2 -,pi,P2) = J exp{- 2 i(pija + P23/2)} 

x \f)*(x x + y 1 , x 2 + y 2 )tf>(xi - y\,x 2 - y 2 )dy\dy 2 . (55) 

The Wigner function corresponding to the oscillator wave function of Eq.(51) is 

/1\ 2 f a . a 2 

W(x l ,x 2 ;pi,p 2 ) = (-J exp |-e T? (ar 1 cos- - x 2 sm -) 

—e~ n (xi sin ^ + x 2 cos ^) 2 - e~ v (p t cos ^ - p 2 sin ^) 2 

-e T '(pi sin ^ + p 2 cos ^) 2 | . (56) 

If we do not make observations in the x 2 p 2 coordinates, the Wigner function becomes 

M / (x,,pi)= W(x u x 2 -,p uP2 )dx 2 dp 2 . (57) 

The evaluation of the integral leads to 

1 


1/2 


W(x uPl ) = 


7r 2 (l + sinh 2 t] sin 2 a) 


£h |_ Pi 

cosh r/ — sin y cos a cosh y + sin y cos a 

This Wigner function gives an elliptic distribution in the phase space of x, and p } . This distribution 
gives the uncertainty product of 

(Ax) 2 (Ap) 2 = i(l + sinh 2 y sin 2 a). (59) 

This expression becomes 1/4 if the oscillator system becomes uncoupled with o = 0. Because Xj 
is coupled with x 2 , our ignorance about the x 2 coordinate, which in this case acts as Feynman s 
rest of the universe, increases the uncertainty in the X] world which, in Feynman s words, is the 
system in which we are interested. 

In the Wigner phase-space picture, the uncertainty is measured in terms of the area in phase 
space where the Wigner function is sufficiently different from zero. According to the Wigner 
function for a thermally excited oscillator state, the temperature and entropy are also determined 
by the degree of the spread of the Wigner function phase space. 


x exp 
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